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(57) Abstract 



A device (10, 30, 50) includes a microprocessor or microcontroller and measures real and imaginary parts of complex immittance of 
a cell or battery (20) at n discrete frequencies, where n is an integer number equal to or greater than 2. The device determines cell/battery 
properties by evaluating components of an equivalent circuit model comprising 2n frequency-independent elements. Equating real and 
imaginary parts of measured immittance (60) to values appropriate to the model at the n measurement frequencies (70) defines a system 
of 2n nonlinear equations. Introducing 2n intermediate variables permits solving these equations and leads to values for the 2n model 
elements. A table of element values contains virtually the same information as the spectrum of complex immitance over a wide frequency 
range but provides this information in a more concise form that is easier to store analyze, and manipulate. Thus, the 2n element values 
may themselves comprise the desired result. 
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METHOD AND APPARATUS FOR DETERMINING 
BATTERY PROPERTIES FROM COMPLEX 
IMP EDANCE / ADMITTANCE 

BACKGROUND OF THE INVENTION 
5 Small -signal ac measurement techniques have 

proven useful in determining properties of 
electrochemical cells and batteries such as cranking 
power, percent capacity, and state-of -health . These 
techniques have generally utilized single- frequency 
10 measurements of a single quantity, such as 
conductance (e.g., U.S. patents 5,585,728 and 
5,140,269 to Champlin) , resistance (e.g., U.S. patent 
3,676,770 to Sharaf et al, U.S. patent 3,753,094 to 
Furuishi, U.S. patent 5,047,722 to Wurst et al) , or 
15 "impedance" (e.g., U.S. patent 4,697,134 to Burkum et 
al, U.S. patent 5,773,978 to Becker). However, 
considerably more information of an electrical, 
chemical, and physical nature is contained in the 
continuous spectrum of complex immittance, i.e., 
20 either impedance or admittance, displayed over a 
range of frequencies. (See, e.g., David Robinson, 
: ^Electrochemical Impedance Spectroscopy in Battery " 
Development and Testing", BATTERIES INTERNATIONAL, 
31, pp. 59-63, April, 1997). A big challenge for 
25 field testing batteries is to acquire such 
information from a relatively small number of 
measurements obtained at a few selected "spot" 
frequencies. 

Muramatsu discloses one -approach to this 
30 challenge in U.S. patent 4,678,998. He measures 
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impedance magnitude at two frequencies. At each 
frequency he compares the measured magnitude with 
that of a predetermined experimental relationship 
between impedance magnitude, remaining capacity, and 
5 remaining service life. He reports that such 

measurements can separately determine the battery' s 
remaining capacity and its . remaining service life. 
Randin discloses a second approach in U.S. Patent 
4,743,855. He reportedly determines a battery's 

10 state-of -discharge from the argument (i.e., phase 
angle) of the difference between complex impedances 
measured at two frequencies. Bounaga discloses still 
another approach in U.S. patent 5,650,937. He 
reportedly determines state-of -charge from 

15 measurements of only the imaginary part of complex 
impedance obtained at a single frequency. All three 
of these approaches have fairly limited objectives, 
however. Much more information is actually contained 
in the complete spectrum of complex immittance than 

20 is acquired by Muramatsu, Randin, or Bounaga. 

Equivalent circuit modeling may assist one 
in relating complex immittance spectra to electrical, 



chemical, or physical properties of a battery. A 
complex nonlinear least -squares curve- fitting 

25 procedure has been used by electrochemists to relate 
impedance spectra to nonlinear electrochemical 
models. (See, e.g., J. Ross Macdonald and Donald R. 
Franceschetti, "Precision of Impedance Spectroscopy 
Estimates of Bulk, Reaction Rate, and Diffusion 

30 Parameters", Journal of Electroana.lytica.1 Chemistry, 
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307, pp. 1-11, 1991; see also Bernard A. Boukamp, "A 
Package for Impedance/Admittance Data Analysis", 
Solid State Ionics, 18, pp. 136-140, 1986). This 
complex procedure, however, requires measuring the 
5 complete spectral distribution of cell/battery 
impedance and then making initial estimates of the 
model's parameters to ensure ultimate convergence. 

An equivalent circuit model is an 
interconnection of electrical elements introduced to. 
10 represent terminal characteristics of the battery. 
In a linear small-signal model, these elements 
comprise discrete resistances capacitances and 
inductances. Such models have been described by a 
number of workers including Hampson, et al (N. A . 
15 Hampson, et al , u The Impedance of Electrical Storage . 
Cells", Journa.1 of Applied Electrochemistry, 10, 
pp. 3-11, 1980), Willihnganz and Rohner (E. 
Willihnganz and Peter Rohner, "Battery Impedance", 
Electrical Engineering, 78, No. 9, pp. 922-925, 
20 September, 1959), and DeBardelaben (S. DeBardelaben, 
"Determining the . End of Battery Life", INTELLEC 86, 
IEEE Publication CH2328-3/86/0000-0365 , pp. 365 - 
T8~6~ 1986; and S~. DeBardelaberT; ""A Look aE the~~™ ™~ 
Impedance of a Cell", INTELLEC 88, IEEE Publication 
25 CH2S53-4/88/000-0394 , pp. 394 - 397, 1988) . However, 
none of these workers has disclosed means for 
determining component values of an equivalent circuit 
model from a small number of measurements obtained at 
a few selected u spot" frequencies. That is an 
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important contribution of the invention disclosed 
herein. 

SUMMARY OF THE INVENTION 
A device includes a microprocessor or 
5 microcontroller and measures real and imaginary parts 
of complex immittance of a cell or battery at n 
discrete frequencies, where n is an integer number 
equal to or greater than 2 . The device determines 
cell/battery properties by evaluating components of 

10 an equivalent circuit model comprising 2n frequency- 
independent linear electrical elements. Equating 
measured real and imaginary parts of complex 
immittance to theoretical real and imaginary values 
appropriate to the model at each of the n discrete 

15 frequencies defines a set of 2n nonlinear equations 
in 2n unknowns. By introducing 2n intermediate 

variables, this formidable problem is made linear and 
is systematically solved for the values of the 
components of the model. Once these values are 

20 known, a table of the 2n element values contains 
virtually the same information as the continuous 
spectrum of complex immittance displayed over a range 
of frequencies . However, the rable of values" 

provides this information in a much more concise form 

25 that is easier to store, analyze, and manipulate. 
Thus, circuit element values may themselves comprise 
the desired result. Moreover, the circuit elements 
represent actual processes occurring within the 
battery. Accordingly, a predetermined relationship 

30 between one or more of the elements and an additional 
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electrical, chemical, or physical property of- the 
cell/battery may be invoked to determine the 
additional property. 

The method and apparatus disclosed herein 
5 are efficient, accurate, and easily implemented with 
^" a microcontroller or microprocessor. The invention 
is suitable for a variety of diagnostic applications 
ranging from hand-held battery testers to "smart" 
battery chargers and battery "fuel gauges" in 
10 electric vehicles. Although a lead-acid automotive 
storage battery is used as an example to teach the 
method, the invention is equally applicable to both 
primary and secondary cells and batteries, and to 
those employed in a variety of other applications 
15 and/or employing other chemical systems . 

BRIEF DESCRIPTION OF THE DRAWINGS 
FIG. la is a spectral plot of measured real 
part of admittance of an actual 12 -volt, automotive 
.storage battery. 
20 FIG. lb is a spectral plot of measured 

imaginary part of admittance of an actual 12 -volt 
automotive storage battery. 

FIG. 2 depicts a general small - signal - ~~ ac 
equivalent circuit model of a cell or battery 
25 comprising 2n frequency- independent linear elements. 

FIG. 3 depicts the equivalent circuit model 
of FIG. 2 with n=2 . 

FIG. 4a is a plot of the experimental data 
of FIG. la compared -with a theoretical curve 
30 calculated from the model of FIG. 3. 
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FIG. 4b is a plot of -the experimental data 
of FIG . lb compared with a theoretical curve 
calculated from the model of FIG. 3. 

FIG. 5 depicts the equivalent circuit model 
5 of FIG. 2 with n=3 . 

FIG. 6a is a plot of the experimental data 
of FIG. la compared with a theoretical curve 
calculated from the model of FIG. 5. 

FIG. 6b is a plot of the experimental data 
10 of FIG. lb compared with a theoretical curve 
calculated from the model of FIG. 5. 

FIG. 7 is a block diagram of a device for 
determining battery properties from spot - frequency- 
complex immittance according to the present 
15 invention. 

FIG. 8 is a flow chart depicting the 
control algorithm for the invention embodiment 
disclosed in FIG. 7. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 
20 The impedance of a cell or battery is a 

complex quantity. At a particular discrete or v spot" 
frequency f k , the complex impedance can be written in 
terms of~^L"ts real and~~imaginary parts as 

25 Z(f k ) = R(f k ) + jX(f k ) (1) 

where j = V— T . The real quantities R(f k ) and X(f k ) are, 
respectively, the resistance and the reactance of the 
cell/battery at the frequency f k . They physically 
30 represent ratios of in -phase voltage amplitude to 
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current amplitude, and quadrature voltage amplitude 

to current amplitude, respectively, at- the frequency 

The admittance or a cell or battery is 
likewise a complex quantity. At a particular 

discrete or "spot" frequency f*, the complex 
admittance can be written 

Y(f k ) = G(f k ) + jB(f k ). (2) 

The real quantities G(f k ) and B(f k ) are, respectively, 
the conductance and the susceptance of the 
cell/battery at the frequency f*. They physically 
represent ratios of in-phase current amplitude to 
voltage amplitude, and quadrature current amplitude 
to voltage amplitude, respectively, at the frequency 
fk- 

Complex admittance and complex impedance 
are related to each other by the reciprocal- 
relationship 



Y(.f lr )^i/z(.f^_. a>-- 

Accordingly, spectra of complex admittance and 
spectra of complex impedance contain exactly the same 
information about the cell or battery. The term 
"immittance" will herein denote either quantity in 
instances where the choice is immaterial. 

Measured spectra of real and imaginary 
parts of complex admittance of ■ a typical automotive 



WO 00/16083 



P C T/U S 9 9/2 08^8 



-8- 

storage battery are disclosed in FIGS, la and lb, 
respectively, over the frequency range from 5 Hz to 
1000 Hz. Considerable information about a battery is 
expressed in such spectral plots. One sees from FIG. 
5 lb, for example, that the battery passes through 
series resonance near 250 Hz, being capacitive (B>0) 
below this frequency and inductive B<0) above. 
However, most of the battery information displayed in 
FIGS. la and lb is very subtle, and not at all 

10 obvious from the plots. 

FIG. 2 discloses a small-signal ac 
equivalent circuit model introduced herein to assist 
in reducing spectral plots of complex immittance, 
such as those displayed in FIGS, la and lb, to a 

15 small set of frequency- independent parameters. One 
sees that the model of FIG. 2 comprises a series 
interconnection of a single two-element series R-L 
subcircuit and n-1 two-element parallel R-C 
subcircuits. For n=2 , this general equivalent 

20 ' circuit model reduces to the simple model discussed 
by both Willihnganz and Rohner and by DeBardelaben . 

A complete disclosure of my method for 
"determining circuit -model element values ""from 
measured values of spot -frequency complex immittance 

25 follows. I begin with an expression for the complex 
impedance of the model of FIG. 2 



1 

Z = R + jX = Rl + jcoLl-f — + + ■ 



1/R2 + jcoC2 1/Rn + jwCn 



(4) 
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where co = 27if is the angular frequency. The 

inductance and the capacitances are eliminated from 
this expression by writing them in terms of time 
constants 



x, =L1/R1 
x 2 = R2C2 



r„ =RnCn 



(5) 



The result is 

Z = R + jX = Rl(l + jcox 1 )-h A R2 + + . (6 ) 

(l + j©xj (1 + jcoxJ 

Multiplying both sides of equation (6) by the product 
(l-hja>x 2 ) (l + jcox n ) clears the fractions and yields 



(R + jX)Xl + j©T 2 )-(l + ^ (?) 
R2(l + jcoT 3 )• • • (l + jcoT n ) + • . • + Rn(l + jcox 2 )• • (l + joox^ ) * 

15 

For n=3, equation (7) reduces to " 



(R + jXXl + " jCOT 2 X 1 + j°> T 3 ) = Rl(l + X 1 + J<° T 2 X 1 + ) ( g 

+ R2(l + jcox 3 )+ R3(l + jcox 2 ) * 



20- Equation (7) is divided into two equations by 
multiplying the terms out, separating them into real 
and imaginary parts, and then separately equating 
real to real, and imaginary to imaginary. For n=3, 
this procedure leads to 
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Real Part : 

((D 2 R)(T 2 T3)-h(c0XXT 2 +X3)-03 2 {Rl(T 2 T3+T3T 1 + T,T 2 )} (g) 

+ (Rl-hR2 + R3)=R 

5 Imaginary Part : 

(co 2 xXx 2 t 3 )-(q)RXt 2 + t 3 )+g>{R1(t i +t 2 + x 3 ) + R2t 3 + R3x 2 } ^ 

-CD J ^lfcT^ )}«X 

Equations (9) and (10) are nonlinear since 
10 the 2n unknown resistances and time constants appear 
as combinations of products. I attack this 

formidable problem by defining a new set of 2n 
intermediate variables. These new variables are the 
various combinations of the model's resistances and 
15 time constants that multiply functions of battery 
resistance, battery reactance, and frequency. For 
n=3, the six intermediate variables are defined by 



=(x 2 +x 3 ) (Ha) 
T 2 ^(t 2 x 3 ) (lib) 



*F 3 =(R1 + R2 + R3) (He) 

=(tj +t 2 +t 3 )R1 + t 3 R2 + x 2 R3 (lid) 

T 5 -(x 2 x 3 +x 3 x 1 +x l x 2 )Rl (lie) 

^ S (W3>M < llf > 

20 

When expressed in terms of these new 

variables, equations (9) and (10) are linear. At the 

angular spot frequency (0 k , these two equations are 
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{co k X(o k )}^ + fc*R((D k )>F 2 4- {l}T 3 + {0}^ 4 - {co k K + Wf 6 = R(co k ) 

(12) 

and 

5 - {co k R(co k )}% + fox(© k ))"F a ■+ {OK ■ + {co k }*F 4 + {o}* 5 - {co 3 k }? 6 = X(co k ) - 

(13) 

In general, equations such as (12) and (13) 
comprise a pair of linear inhomogeneous equations for 

1G the 2n intermediate variables, *F, , '-^^n • Even 
though linear, such equations are still not solvable 
since they contain an insufficient number of 
relationships between variables. However, by 

evaluating the complex impedance at n discrete 

15 frequencies, C0j co n , the two equations expand into a 

solvable set of 2n linear inhomogeneous equations in 
2n unknowns . ' Such a system can be solved by the 
well-known method known as Cramer's rule. Cramer's 

rule expresses the 2n solutions, ,\P 2n , as ratios 

20 of determinants having 2n columns and 2n rows. For 
n=3, these six solutions are of the .form 



V^AjAv; ;*F 6 =A 6 /A D (14) 

25 where A D and A,---A 6 *are (6 x 6) determinants given by 



T 
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(IS) 



(2 0) 



(21) 



The determinants disclosed in equations 



(15) - (21), can be systematically evaluated from spot- 
10 frequency immittance measurements by well-known 
numerical expansion techniques. Once their values 

are known, the intermediate variables 4*,, , V F 6 follow 

from equations (14) . ■ The defining equations of the 
intermediate variables, equations (11a) - (llf ) , are 
15 then combined in particular ways to evaluate the 
elements of the equivalent circuit model . 
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The technique proceeds as follows. One 
first combines the defining equations for n-1 of the 
2n intermediate variables to yield an equation for 

the n-1 capacitive time constants T 2 ,---,T n . These n-1 
5 equations are identified by not containing 
resistances. For example, for n=3, I combine 

equation (11a) and equation (lib) to obtain the 
following quadratic equation: 



10 < 3 -^,3+^2=0. (22) 

The two roots of equation (22) are given by the well- 
known quadratic formula 



15 



2 , 3 = ^±^/2) 2 -% - (23) 



For the general case of arbitrary n, combining the n- 
1 defining equations that contain no resistances 
leads to the following polynomial equation of order 
20 n-1: 



x (- 1 ).y iX (-2) + ± xj/ ni =0 (24) 

The n-1 roots of equation (24) are the capacitive 
25 time constants t 2 ,---,T n . Although general formulas 
similar to equation (23) do not exist to solve higher 
order polynomial equations, the roots of equation 
(2 4) can always be found using well-known numerical 
root - finding algorithms - 
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Once the capacitive time constants have 
been determined, the inductive time constant t, 
follows by eliminating Rl from the two defining 
equations for intermediate variables that are 

5 proportional to Rl . For example, with n=3 , I combine 

equations (lie) and (llf) to obtain 

T,={K/^ 6 )-iA 2 -i/T 3 r. (25) 

10 For the general case, of arbitrary n, the expression 
for x 1 is of the form 

*l ={K„-,)/ l * / 2n )-lA 2 - -lAn}"' • (26) 

15 Thus, in principle, all n time constants are known. 

In general, there are 2n equations defining 
the 2n intermediate variables- The first n-1 of 
these equations contain no resistances and are 
employed to determine the capacitive time constants. 

20 Two of the remaining n+1 equations are employed to 
determine the inductive time constant . By choosing 
either one- of these two equations, along with all of 
the remaining n-1 unused equations, one obtains a set 
of n linear inhomogeneous equations in the n unknowns 

25 Rl,- -,Rn . Cramer's rule can then be invoked to solve 
this system for the values of these n resistances. 
For example, with n=3 / I use equations (16c), (16d) , 
and (16f) to obtain the following set of three linear 
equations in three unknowns: 
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{1}R1 

{tj + x 2 +t 3 }R1 



x 2 x 3 }Rl 



+ {1}R2 
{x 3 }R2 
+ {0}R2 



+ {l}R3 = ¥ 3 
+ {t 2 }R3 = ¥ 4 
+ {0}R3 = ¥ 6 



(27) 



Cramer's rule yields the following three solutions: 



Rl = 





^3 


1 


1 






*4 


^3 








% 


0 


0 






1 




1 


1 


x, + x 2 +x 3 


T 3 




x,x 2 x 3 




0 


0 



(28) 



and 



R2 = 



x, + x 2 +x 3 



^4 



1 



x, + x 2 +x. 



(29) 





1 1 % 

x,+x, +x 3 x 3 % 




R3 = 


x,x,x 3 0 T 6 

1 1 1 

x,+x, +x 3 x 3 - x 2 
XjX 2 x 3 0 0 


(30) 



Thus, in principle, all n resistance values are 
known . 



10 
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Finally, I invert equations (5) to 
determine the single inductance value and the n-1 
capacitance values from the n known time constants 
and the n known resistance values 



LI = t,R1 (31a) 
C2 = x 2 /R2 (31b) 

Cn = i n /Rn (31c) 

•This completes the disclosure of. my method 
for determining the 2n element values. Although n is 
unrestricted in principle, the problem of expanding 
large determinants will probably limit n to the range 
2<n<8. In summary, the steps are: 

1. One first finds 2n intermediate 
variables by evaluating 2n ratios of 
15 determinants comprising 2n columns and 2n rows. 

The determinants comprise sums and differences 
of products combining the n spot frequencies 
with real and imaginary parts cf impedance or 
admittance at the n spot frequencies. 
"20 ; 2~ The n^l capacifive ETme constants are" " 

found as roots of a polynomial equation of order 
n-1. The polynomial's coefficients comprise n-1 
of the 2n intermediate variables whose defining 
equations contain no resistances. 
25 3. The inductive time constant is 

determined . from an equation containing the 
capacitive time constants found in step 2 along 
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with the ratio of the two intermediate variables 
proportional to Rl . 

4. The n resistances are determined by 
evaluating n ratios of (n x n) determinants. 

5 Elements of these determinants contain the n 

time constants found in steps 2 and 3 along with 
n of the remaining n+1 intermediate variables 
not employed in step 2 . 

5. Finally, the one inductance and the n- 
10 1 capacitances are evaluated by combining each 

of the ii resistances with a corresponding time 
constant . 

In practice, this procedure is readily implemented in 
software . 

15 The exact procedure disclosed above can be 

simplified by choosing one of the n spot frequencies, 

say Q) n , to be sufficiently high that the impedance of 
the series chain of R-C subcircuits is negligibly 
small at this frequency. A preliminary measurement 

20 of Z(co n ) then gives approximations to RI and LI 

directly. Subtracting Rl + jco k Ll from the total 

mea sured impedance Z( co k ) at each of the n-1 remaining 

spot frequencies and equating this result to the 
theoretical impedance of the R-C subcircuit chain at 

25 each frequency leads to a system of 2n-2 equations in 
2n-2 unknowns. Solving for the appropriate 2n-2 
intermediate variables by Cramer's rule involves 
evaluating determinants that are fewer in number and 
smaller in size than the determinants evaluated in 

30 the exact procedure disclosed above. 
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FIG. 3 discloses an equivalent circuit 
model derived from FIG. 2 by letting n=2 . This 
simple model is essentially the one discussed by 
Willihnganz and Rohner and by DeBardelaben in the 

5 publications cited earlier. By using the exact _ 
procedure" disclosed above, I evaluated the four 
linear circuit elements of the model of FIG. 3 from 
the experimental data disclosed in FIGS, la and lb at 
the two spot frequencies £ a = 5 Hz and f 2 =. 1000 Hz. 

10 The results .of this evaluation procedure are 
displayed in Table 1. 



15 



20 



25 



Table 1. 



Model Element Values for n=2 



30 



Rl = 4.388 mQ R2 = 12.987 mQ 

LI = 0.3885 HH C2 - 2.602 F 

Theoretical curves of the real and 
imaginary parts of admittance as functions of 
frequency were calculated for the model of FIG. 3 by 
assuming the element values displayed in Table 1. 
The theoretical curves are plotted along with the 
measured curves for comparison in FIGS 4a and 4b. 
One sees that the experimental and theoreticarcurveT 
agree exactly at the two spot frequencies as would be 
anticipated. However, away from the spot frequencies 
the agreement is seen to be poor. This indicates that 
the model of FIG. 3 does not adequately represent the 
battery over this frequency range. 

Much better agreement is obtained with the 
n=3 model depicted in FIG. 5. Using the exact 
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procedures disclosed above, I evaluated the six 
linear circuit elements of the model of FIG. 5 from 
the experimental data disclosed in FIGS, la and lb at 
the three spot frequencies fi = 5 Hz, f 2 = 7 0 Hz, and 
5 f 3 = 1000 Hz. The results are displayed in Table 2. 



Table 2. Model Element Values for n=3 
Rl = 4.381 mD R2 = 1.227 mQ R3 = 13.257 jdQ 
10 LI = 0.411 JJ.H C2 = 1.812 F C3 = 3.14 F 

Theoretical curves of the real and 
imaginary parts of admittance as functions of 
frequency were calculated for the model of FIG. 5 by 
15 assuming the element values displayed in Table 2 . 
These curves are plotted along with the measured 
curves for comparison in FIGS. 6a and Sh . Once again 
one sees that the experimental and theoretical curves 
agree exactly at the spot frequencies. Away from 

20 these frequencies, however, the agreement is seen to 
also be very good. Such good agreement proves the 
model of FIG. 5 to be an excellent representation of 
the battery over the frequency range from 5 Hz to 
1000 Hz. Accordingly, the n=3 model much more 

25 closely describes actual processes within the battery 
than does the n=2 model. 

The excellent agreement between the 
experimental curves and the theoretical predictions 
of the model means that Table 2 contains virtually 

30 the same information about' the battery as does the 
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continuous spectrum of complex admittance displayed 
in FIGS, la and lb. However, Table 2 provides this 
information in a much more concise form that is far 
easier to store, analyze, and manipulate. 
5 Accordingly, the information displayed in Table 2 may 

Ttself comprise ' CEeT desired result:". 

Moreover, since the circuit elements 
defined in the extended model closely describe actual 
processes occurring within the battery, a 
10 predetermined relationship can additionally be 
invoked if desired, • to implement a final step of 
determining one or more additional battery property. 
For example, I have found that the battery's cold- 
cranking ampere (CCA) capacity is quite accurately 
15 given by 



CCA = 2662/R1 (32) 

where Rl is expressed in milliohms. Thus, the 

20 battery whose complex admittance spectrum is 
disclosed in FIGS, la and lb is capable of supplying 
608 cold-cranking amperes. This important 

information is not at all obvious from the spectral 
plots of FIGS. la and lb. Other electrical 

25 properties such as state-of -charge and ampere-hour • 
capacity; chemical properties such as ion 
concentration and plate composition; and physical 
properties such as battery temperature and effective 
plate area, find similar expression in the complex 
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immittance spectra of the battery and may be 
determined in a comparable manner. 

Figure 7 discloses a block diagram of a 
device for determining cell/battery properties from 
5 spot -frequency complex immittance according to the 
present invention. Measuring circuitry 10 

electrically couples to cell/battery 2 0 by means of 
current -carrying contacts A and B and voltage- sensing 
contacts C and D. Measuring circuitry 10 passes a 

10 periodic time-varying current i(t) through contacts A 
and B and senses a periodic time -varying voltage v(t) 
across contacts C and D. By appropriately processing 
and combining i(t) and v(t) , measuring circuitry 10 
determines real and imaginary parts of complex 

15 immittance at a measuring frequency f k ; where f k is a 
discrete frequency component of waveforms i(t) and v(t) . 

Control circuitry 3 0 couples to measuring 
circuitry 10 via command path 4 0 and commands 
measuring circuitry 10 to determine the complex 

20 immittance of cell/battery 20 at each one of n 
discrete measuring frequencies, where n is an integer 

n umber equal to or great er than 2. This action 

defines 3n experimental quantities: the values of the 
n measuring frequencies and the values of the n 

25 imaginary and n real parts of the complex immittance 
at the n measuring frequencies. 

Computation circuitry 50 couples to 
measuring circuitry 10 and to control circuitry 3 0 
via data paths 60 and 70, respectively, and accepts 

30 . the 2n experimental values from measuring circuitry 
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10 and the values of the n measuring frequencies from 
control circuitry 30. Upon a u Begin Computation" 
command from control circuitry 3 0 via command path 
80, computation circuitry 50 invokes the procedure 
5 disclosed above to combine these 3n quantities 
numerically / to evaluate the 2n elements of the 
equivalent circuit model. These 2n element values 
may themselves comprise the desired result- However, 
if "desired, computation circuitry 50 can also perform 
10 an additional step by relating one or more of the 
model element values to an additional cell/battery 
property to determine the additional property. 

In "practice, a single microprocessor or 
microcontroller running an appropriate software 
15 program can perform the functions of both control 
circuitry 3 0 and computation circuitry 50 as well as 
much of the function of measuring circuitry 10. 
Microprocessor controlled impedance measuring 
apparatus is disclosed in a copending U.S. patent 
20 application. 

Figure 8 discloses a flow chart depicting a 
control algorithm for the invention embodiment of 
FIG. 7. Upon entering the procedure at 100, control 
circuitry .30 initializes a counter at 105 used to 
25 identify each spot frequency .f k . At 110, control 
circuitry 3 0 commands measuring circuitry 10 to 
excite the cell/battery with a periodic signal having 
a sinusoidal component at frequency f k and to 
determine the real and imaginary parts of complex 
30 immittance at this frequency. At decision block 115, 
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control circuitry 3 0 determines whether this 
procedure has been performed at all of the desired 
spot frequencies. If the answer is no, the counter 
is incremented at 12 0 and the procedure repeated at a 
5 new frequency. If yes, control circuitry 3 0 commands 
computation circuitry 50 to begin the computation at 
125. Computation circuitry 50 begins at step 130 by 
determining the 2n intermediate variables from the 
values of real and imaginary parts of complex 

10 immittance at the n spot frequencies along with the n 
spot frequency values themselves. The 2n frequency- 
independent model elements are then evaluated from 
the 2n intermediate variable values at step 135. 
Finally, at optional step 14 0 computation circuitry 

15 50 can, if desired, invoke an additional 
predetermined relationship between one or more of the 
model elements and a desired cell/battery property to 
determine the desired property. 

This completes the disclosure of my 

20 invention. The method and apparatus are efficient, 
accurate, and easily implemented with a 

microcontroller or microprocessor. The invention is 
quite general and stTitable for — a — wdrde — variety — of- 
diagnostic applications ranging from hand-held 

25 battery test instruments to "smart" battery chargers 
and battery n f uel gauges" in electric vehicles . 
Although a lead-acid storage battery was used as an 
example to teach the method, the disclosed invention 
is equally applicable to both primary and secondary 

30 cells and batteries, and to cells/batteries employed 
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in a variety of other applications and/or employing 
other chemical systems . 

The present invention has been described 
with reference to a preferred embodiment. However, 
5 workers skilled in the art will recognize that 
changes may be made in form and detail without 
departing from the spirit and scope of the invention. 
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WHAT IS CLAIMED IS: 

1. A device for determining values of at least 

one of 2n elements comprising an equivalent circuit 
model of an electrochemical cell or battery where n 
5 is an integer equal to or greater than two, said 
device comprising: 

measuring circuitry adapted to couple to 
said cell or battery and adapted to 
pass a periodic current through said 
10 cell or battery, to sense a periodic 

voltage across said cell or battery, 
and to determine a real part and an 
imaginary part of complex immittance 
of said cell or battery at a 
15 measurement frequency comprising a 

component frequency of said periodic 
current and said periodic voltage ; 
control circuitry coupled to said 
measuring circuitry and adapted to 
20 command said measuring circuitry to 

select each one of n said measurement 
frequencies thereby defining n said 
real - parts"; n~ sa"i _ d-d:ma-g±na-ry--pa-rtrs— -and— 
n said measurement frequencies; and, 
25 computation circuitry coupled to said 

measuring circuitry and to said 
control circuitry and adapted to 
numerically combine values of said n 
real parts, said n imaginary parts, 
30 and said n measurement frequencies to 
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determine said values of at least one 
of said 2n elements comprising said 
equivalent circuit model of said 
electrochemical cell or battery. 

5 

2 i A device as in claim 1 wherein n is equal 
to three . 

3 . A device as in claim 1 wherein said control 
10 circuitry and said computation circuitry comprise a 

mioroprocessor or microcontroller running a software 
program adapted to select each one of said n 
measurement frequencies and to numerically combine 
values of said n real parts, said n imaginary parts, 
15 and said n measurement frequencies to determine said 
values of at least one of said 2n elements comprising 
said equivalent circuit model of said electrochemical 
cell or battery. 

20 4 . A device as in claim 3 wherein said 

software program is further adapted to numerically 
combine values of said n real parts, said n imaginary 
parts, and said n measurement frequencies to evaluate" 
2n intermediate variables, and said values at least 

25 one of said 2n elements are determined from values of 
said 2n intermediate variables . 

5. A device for determining an electrical, 

chemical, or physical property of an electrochemical 
30 cell or battery comprising: 
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measuring circuitry adapted to couple to 
said cell or battery and adapted to 
pass a periodic current through said 
cell or battery, to sense a periodic 
5 voltage across said cell or battery, 

and to determine a real part and an 
imaginary part of complex immittance 
of said cell or battery at a 
measurement frequency comprising a 
10 component frequency of said periodic 

current and said periodic voltage; 
control circuitry coupled to said 
measuring circuitry and adapted to 
command said measuring circuitry to 
15 select each one of n said measurement 

frequencies thereby defining n said 
real parts, n said imaginary parts and 
n said measurement frequencies where n 
is an integer number equal to or 
20 greater than two; and, 

computation circuitry coupled to said 
measuring circuitry and to said 
— - control circuitry and "adapted to 



25 



numerically combine values of said n 
real parts, said n imaginary parts, 
and said n measurement frequencies to 
determine said electrical, chemical, 
or physical property of said 
electrochemical cell or battery. 



30 
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6. A device as in 'claim 5 wherein n is equal 
to three . 

7. A device as in claim 5 wherein said control 

5 circuitry and said computation circuitry comprise a__ 
microprocessor or microcontroller running a software 
program adapted to select • each one of said n 
measurement frequencies and to numerically combine 
values of said n real parts, said n imaginary parts, 
10 and said ■ n measurement- frequencies to determine said 
electrical, chemical, or physical property of said 
electrochemical cell or battery. 

8. A device as in claim 7 wherein said 
15 software program is further adapted to numerically 

combine values of said n real parts, said n imaginary 
parts, and said n measurement frequencies to evaluate 
2n intermediate variables, and said electrical, 
chemical, or physical property of said 

20 electrochemical cell or battery is determined from 
values of said 2n intermediate variables . 



9. a device as in claim 8 wherein said 

software program is further adapted to numerically 
25 combine said 2n intermediate variables to evaluate 
circuit model elements and said electrical, chemical, 
or physical property of said electrochemical cell or 
battery is determined from values of said circuit 
model elements . 
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10. A method for determining an electrical, 

chemical, or physical property of an electrochemical 
cell or battery comprising the steps of: 

measuring real and imaginary parts of 
5 complex immittance of said 

electrochemical cell or battery at n 
discrete frequencies where n is an 
integer number equal to or greater 
than two ; 

10 evaluating 2n intermediate variables by 

numerically combining values of said n 
discrete frequencies and values of 
said real and imaginary parts of said 
complex immittance at said n discrete 
15 frequencies ; 

evaluating n-1 capacitive time constants by 
numerically combining values of n-1 of 
said 2n intermediate variables; 
evaluating an inductive time constant by 
20 numerically combining values of said 

n-1 capacitive time constants and 
values of two of said 2n intermediate 
~ — — variables ; 



evaluating n resistive elements by 
25 numerically combining values of said 

n-1 capacitive time constants, the 
value of said inductive time constant, 
and values of n of said 2n 
intermediate variables; 
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evaluating n reactive elements by 
numerically combining each value of 
said n resistive elements with a 
corresponding value of one of said n-1 
capacitive time constants and said 
inductive time constant; and 
relating values of one or more of said 
resistive elements and said reactive 
elements to said electrical, chemical, 
or physical property of said 
electrochemical cell or battery. 

11- A method as in claim 10 wherein n is equal 

to three. 



12 • A method for determining element values of 

an equivalent circuit model of an electrochemical 
cell or battery, said model comprising a series 
interconnection of a single two-element series R-L 
20 subcircuit and n-1 two-element parallel R-C 
subcircuits where n is an integer number equal to or 
greater than two, said method comprising the steps 



Of: 

measuring real and imaginary parts, of 
25 complex immittance of said 

electrochemical cell or battery at n 
discrete frequencies; 
evaluating 2n . intermediate variables by 
numerically combining values of said n 
30 discrete frequencies and values of 
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said real and imaginary parts of said 
complex immittance at said n discrete 
frequencies ; 
evaluating n-1 capacitive time constants by 
numerically combining values of n-1 of 
said 2n intermediate variables ; 
evaluating an inductive time constant by 
numerically combining values of said 
n-1 capacitive time constants and 
values of two of said 2n intermediate 
variables ; 

evaluating a resistance of said single R-L 
subcircuit and n-1 resistances of said 
n-1 R-C subcircuits by numerically 

15 combining values of said n-1 

capacitive time constants, the value 
of said inductive time constant, and 
values of n of said 2n intermediate 
variables; and, 

20 evaluating an inductance of said R-L 

subcircuit and n-1 capacitances of 
said n-1 R-C subcircuits by 



numerically combining a resistance 
value for each said subcircuit with a 
25 corresponding time constant value 

associated with the same said 
subcircuit . 



13. A method as in claim 12 wherein n is equal 

30 to three. 
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14 . A method for determining element values of 

an equivalent circuit model of an electrochemical 
cell or battery, said model comprising .a series 
5 interconnection of a single two-element series R-L 
subcircuit and n-1 two-element parallel R-C 
subcircuits where n is an integer equal to or greater 
than two, said method comprising the steps of: 

measuring real and imaginary parts of 
10 complex immittance of said 

electrochemical cell or battery at n 
discrete frequencies where one of said 
n discrete frequencies is a high 
frequency; 

15 evaluating a resistance and an inductance 

of said R-L subcircuit by numerically 
combining the value of said high 
frequency and values of said real and 
imaginary parts of said complex 

20 immittance at said high frequency; 

evaluating a difference impedance at each 
of n-1 remaining frequencies by 



numerically combining values of said 
resistance, said inductance, and said 
25 remaining frequency, with the value of 

said complex immittance at each said 
remaining frequency; 
evaluating 2n-2 intermediate variables by 
numerically combining values of said 
n-1 remaining frequencies and values 
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of real and imaginary parts of said 
difference impedance at said n-1 
remaining frequencies ; 
evaluating n-1 capacitive time constants of 
5 said n-1 R-C subcircuits by 

numerically combining values of n-1 of 
said 2n-2 intermediate variables; 
evaluating n-1 resistances of said n-1 R-C 
subcircuits by numerically combining 
10 values of said n-1 capacitive time 

constants and values of n-1 of said 
2n-2 intermediate variables; and, 
evaluating n-1 capacitances of said n-1 R-C 
subcircuits by. numerically combining a 
15 resistance value of each said 

subcircuit with a corresponding 
capacitive time constant value of the 
same said subcircuit . 

20 15. A method as in claim 14 wherein n is equal 

to three. 

~~16~i — ^aT" method "for determining el^tricaTlT" 
chemical, or physical properties of an 

25 electrochemical cell or battery comprising the steps 
of : 

measuring real and imaginary parts of 
complex immittance of said 

electrochemical cell or battery at n 
30 discrete frequencies where n is an 
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integer number equal to or greater 
than two and one of said discrete 
frequencies is a high frequency; 
evaluating a first resistance and an 
5 inductance by numerically combining 

the" value" of" said high frequency and 
values of said real and imaginary 
parts of said complex immittance at 
said high frequency; 

10 evaluating a difference impedance at each 

of n-1 remaining frequencies by 
numerically combining values of said 
first resistance, said inductance, and 
said remaining frequency with the 

15 value of ' said complex immittance at 

each of said n-1 remaining 

frequencies ; 
evaluating 2n-2 intermediate variables by 
numerically combining values of said t 

20 n-1 remaining frequencies and values 

of real and imaginary parts of said 
difference impedance at said n-1 



25 



30 



remaining frequencies; 

evaluating n-1 time constants by 
numerically combining values of n-1 of 
said 2n-2 intermediate variables ; 

evaluating n-1 . second resistances by 
numerically combining values of said 
n-1 time constants and values of n-1 
of said 2n-2 intermediate variables; 
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evaluating n-1 capacitances by numerically 
combining the value of each of said n- 
1 second resistances with a 
corresponding value of each of said 
n-1 time constants; and, 

relating values of one or more of said 
first resistance, said n-1 second 
resistances, said inductance, and said 
n-1 capacitances to said electrical, 
chemical, or physical property of said 
electrochemical cell or battery. 

17. A method as in claim 16 wherein n is equal 

to three . 

18 . A device for determining an electrical 

chemical, or physical property of an electrochemical 
cell or battery adapted to perform the steps of the 
method of claim 10. 

19. A device for determining element values of 

an equivalent circuit model for an electrochemical cell 



- — battery, said model comprising a series 
interconnection of single two-element series R-L 
25 subcircuit and n-1 two-element parallel R-G subcircuits 
where n is an* integer number equal to or greater than 
two, said device adapted to perform the steps of: 

measuring real and imaginary parts of 
complex immittance of said 
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electrochemical cell or battery at n 
discrete frequencies ; 
evaluating 2n intermediate variable's by 
numerically combining values of said n 
5 discrete frequencies and values of 

said real and imaginary parts of said 
complex immittance at said n discrete 
frequencies; 
evaluating n-1 capacitive time constants by 
10 numerically combining values of n-1 of 

said 2n intermediate- variables ; 
evaluating an inductive time constant by 
numerically combining values of said 
n-1 capacitive time constants and 
15 values of two of said 2n intermediate 

variables ; 

evaluating a resistance of, said single R-L 
subcircuit and n-1 resistances of said 
n-1 R-C subcircuits by numerically. 
20 combining values of said n-1 

capacitive time constants, the value 
of said inductive time constant, and 
values of n of said 2n intermediate 
variables; and, 
evaluating an inductance of said R-L 
subcircuit and n-1 capacitances of 
said n-1 R-C subcircuits by 
numerically combining a resistance 
value for each said subcircuit with a 
corresponding time constant value 



25 
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associated with the same said 
subcircuit. 

20. A device for determining element values of 

5 an equivalent circuit model of an electrochemical 
cell or battery, said model comprising a series 
interconnection of single two-element series R-L 
subcircuit and n-1 two-element parallel R-C subcircuits 
where n is an integer number equal to or greater than 
10 two, said device adapted to perform the steps of: 

measuring real and imaginary parts of 
complex immittance of said 

electrochemical cell or battery at n 
discrete frequencies where one of said 
15 n discrete frequencies is a high 

frequency; 

evaluating a resistance and an inductance 
of said R-L subcircuit by numerically 
combining the value of said high 
20 frequency and values of said real and 

imaginary parts of said complex 
immittance at said high frequency; 



evaluating a difference impedance at each 
of n-1 remaining frequencies by 
25 numerically combining values of said 

resistance, said inductance, and said 
remaining frequency, with the value of 
said complex immittance at each said 
remaining frequency ; 



WO 00/16083 



PCT/US99/2083S 



-39- 

evaluating 2n-2 intermediate variables by 
numerically combining values of said 
n-1 remaining frequencies and values 
of real and imaginary parts of said 
5 ' difference impedance at said n-1 

remaining frequencies ; 

evaluating n-1 capacitive time constants of 
said n-1 R-C subcircuits by 
numerically combining values of n-1 of 
10 said 2n-2 intermediate variables; 

evaluating n-1 resistances of said n-1 R-C 
subcircuits by numerically combining 
values of said n-1 capacitive time 
constants and values of n-1 of said 
15 2n-2 intermediate variables; and, 

evaluating n-1 capacitances of said n-1 R-C 
subcircuits by numerically combining a 
resistance value of each said 
subcircuit with a corresponding 
20 capacitive time constant .value of the 

same said subcircuit . 



21. A device for determining an electrical, 

chemical, or physical property of an electrochemical 
25 cell or battery adapted to perform the steps of : ' . 

measuring real and imaginary parts of 
complex immittance of said 

electrochemical cell or battery at n 
discrete frequencies where n is an 
30 integer number equal to . or greater 
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than two and one of said discrete 
frequencies is a high frequency; 
evaluating a first resistance and an 
inductance by numerically combining 
5 the value of said high frequency and 

values of said real and imaginary 
parts of said complex immittance at 
said high frequency; 
evaluating a difference impedance at each 
10 of n-1 remaining frequencies by 

numerically combining values of said 
first resistance, said inductance, and 
said remaining frequency with the 
value of said complex immittance at 
15 each of said n-1 remaining 

frequencies ; 
evaluating 2n-2 intermediate variables by 
numerically combining values of said 
n-1 remaining frequencies and values 
20 of real and imaginary parts of said 

difference impedance at said n-1 
remaining frequencies ; 
evaluating n^T" " time - " constants ;By 
numerically combining values of n-1 of 
25 said 2n-2 intermediate variables; 

evaluating n-1 second resistances by 
numerically combining values of said 
n-1 time constants and values of n-1 
of said 2n-2 intermediate variables; 
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-41- 

evaluating n-1 capacitances by numerically 
combining the value of each of said n- 
1 second resistances with a 
corresponding value of each of said 
n-1 time constants; and, 

relating values of one or more of said 
first resistance, said n-1 second 
resistances, said inductance, and said 
n-1 capacitances to said electrical, 
chemical, or physical property of said 
electrochemical cell or battery. 
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